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ABSTRACT 

Quasar feedback outflows are commonly invoked to drive gas out of galaxies in the 
early gas-rich epoch to terminate growth of galaxies. Here we present simulations that 
show that AGN feedback may drive not only gas but also stars out of their host galaxies 
under certain conditions. The mechanics of this process is as following: (1) AGN-driven 
outflows accelerate and compress gas filling the host galaxy; (2) the accelerated dense 
shells become gravitationally unstable and form stars on radial trajectories. For the 
spherically symmetric intial conditions explored here, the black hole needs to exceed 
the host's M a mass by a factor of a few to accelerate the shells and the new stars 
to escape velocities. We discuss potential implications of these effects for the host 
galaxies: (i) radial mixing of bulge stars with the rest of the host; (ii) contribution of 
quasar outflows to galactic fountains as sources of high-velocity clouds; (iii) wholesale 
ejection of hyper velocity stars out of their hosts, giving rise to type II supernavae 
on galactic outskirts, and contributing to reionisation and metal enrichment of the 
Universe; (iv) bulge erosion and even complete destruction in extreme cases resulting 
in overweight or bulgeless SMBHs. 

Key words: galaxies: evolution - quasars: general - accretion, accretion discs - stars: 
hypervelocity stars 



1 INTRODUCTION 

Super-massive black holes (SMBH) accreting gas at rates 
comparable to the Eddington accretion rate are believed to 
power quasars. They are expected to launch powerfu l winds 
i|Shakura fc Sunvaevl [l973l : iKind [20031 ; iProgal 12001 ). Such 
winds are consistent with the fast v ou t ~ 0-lc winds de- 
tected via absorption (e.g., Pounds et a l. 2003a. b]) and also 



recently in emission dPounds fc Vaughan 



I2OIII ) in AGN X- 
ray spectra. These nuclear (most likely R £ 1 pc) winds 
must be wide-angle t o explain their detection frequency 
l|Tombesi et al.l 2010a, b), and energetically should be ca- 
pable of driving outflows clearing significant fracti o ns of 
all g a s of the parent galaxy out l|Silk fc Reesl Il998l ; iKina 
120051 : IZubovas fc Kmefboiiaj h Recent observations detect 
such kpc-scale neutral and ionized outflows with outflow ve- 
locities of ~ 1000 km s _1 and mass outflow rates of hun- 
dreds to thousands of M fl yr -1 (e.g.. iFeruglio et ail l2010l : 
ISturm et ai1l201ll : iRupke fc Veilleuxl |201 ll ) , which are best 
interpreted as the mass-loaded nuclear outflows. These out- 
flows have wi de-ranging implications for galaxy formation 
and evolution jDi Matteo et alj|2003 ). 

Recently, Nayakshi n fc Zubovas! (|2012l , hereafter NZ12) 
have shown that quasar outflows affect the ambient gas in 
the bulge of the host galaxy in two different ways depending 



on whether the ambient shocked gas is able to cool rapidly 
or not. In the gas-poor epoch, e.g., the present day bulge of 
the Milky Way, the shocked gas cools slowly, so the shock 
is approximately adiabatic. The layer of the shocked gas is 
then geometrically thick and the density in it is ~ 4 times 
higher than the pre-shock density, as in the strong-shock 
limit. In contrast, in the gas-rich epoch, when the pre-shock 
gas density in the bulge is higher, the shocked gas cools more 
rapidly than the shock can propagate through the bulge. The 
shocked gas is then swept up into a geometrically thin shell, 
and the density in the shell is roughly v 2 /c 2 times the pre- 
shock density, where v e ~ 1000 km s _1 and c s ^, 10 km s 
are the shock front velocity and the sound speed in the shell, 
respectively. If the pre-shock gas density is a ^ 0.1 fraction 
of the total density responsible for the potential, the swept 
up shell density is orders of magnitude higher than the tidal 
density. Unsurprisingly then, NZ12 found that the shell can 
fragment into star-forming clumps. 

Such positive AGN feedback has been explored be- 
fore with re gard to jet interac t ion with clump y interstel- 
lar medium jGaibler et alj|2012t Isiik et al.ll2012T ). These au- 
thors found that AGN jets can also enhance star formation 
in their host galaxies. However, the details of that process 
are markedly different from what consider here and so we do 
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not attempt to draw direct comparison with these previous 
studies. 

NZ12 concentrated on conditions when the shell stalls 
and falls back on to the SMBH, which phys ically corre s pond s 
to black hole masses below the M a mass (King 2003, 2005) 

M a ~ 3.67 x 10* Moo M S) (1) 

where 0200 is the velocity dispersion in the host galaxy bulge, 
in units of 200 km/s, and / is the ratio of the gas density in 
the bulge to the total density in the bulge, scaled to the cos- 
mological ratio of 0.16. Here we probe the evolution of the 
outflow for Mhh ~ M a - We find that the instability leading 
to the shell fragmentation is still effective when the shell 
is being driven outward at velocities as large as £ 1000 
km s , leading to formation of self-gravitating dense gas 
clumps, and eventually newly born stars, moving radially 
outward with these high velocities. This result suggests that 
quasar outflows could not only be expelling gas from the host 
galaxies but also making some stars by strong compression 
of that material; since the material is streaming away at high 
velocity the new stars are born with that velocity as well. 

Below we present numerical simulations that exemplify 
these effects and also consider their theoretical implications. 



2 NUMERICAL METHODS 

The simulations presented here to demonstrate our main 
point are set up in the manner identical to those presented 
in NZ12, so we only briefly overview the setup. For simplic- 
ity, the host galaxy is modelled by a fixed singular isother- 
mal sphere potential softened in the central 10 pc where 
the SMBH is located. The ambient gas is initialised as a 
sphere with inner radius Ri n — 200 pc and outer radius 
■Rout = lOiiin, at rest, and with the density profile given by 
p g (7?) = fgpo(R) = f s a 2 /2-kGR 2 ; a = 141 km s" 1 is the ve- 
locity dispersion of the potential and f s = 0.16/ is the ratio 
of the gas density to the underlying density of stars and dark 
matter; we choose / = 1 for the simulations presented in this 
paper. We use 5 x 10 5 SPH particles in each simulation, giv- 
ing a mass resolution M TCB ~ 40mspH — 2 x 10 5 M , much 
better than in a typical galaxy-wi de AGN feedback sim ula- 
tions (e.g. Af rcs ~ 6 x 10 6 M in lSpringel et al.ll2005l . and 
much larger than this in cosmological simulations). Initial 
gas temperature is set such that the sound speed is equal 
to a. At time t = 0, the quasar outflow is turned on. The 
critical M a mass above which even a purely momentum- 
driven SMBH outflow should e xpel the gas for this poten- 
tial is M a « 1.5 x 1O 8 /M (jNavakshin fc Powerl l201oh . 
Quasar feedback is implemented via "virtual particles" 
emitted isotropically by the SMBH jN avakshin ct al. 120091 : 
iNavakshin fc Powerllioiol : IZubovas fc Navakshinll2012h . The 
quasar momentum outflow rate is fixed at L^dd/c and the 
energy input into the gas is = 0.05Z/Edd; we use a smooth 
transition between the momentum - and energy-driven flow 
regimes, with the cooling radius l|Zubovas fc King] l2012bl ) 
R c = 500 pc. The SMBH mass is fixed for the duration of a 
simulation; this is a reasonable assumption, given that the 
timescales we are interested in are shorter than the Salpeter 
time (w 40 Myrs). Physically this situation could be realised 
by a small scale dense gas disc feeding the SMBH at close 
to its Eddington rate (small scale gas discs are much more 



resilient to quasar feedback than e xtended spherical distri- 
butions; see INavakshin et al.ll2012h . 

Since we are modelling the inner parts of a single galaxy, 
we can afford a higher numerical resolution, as already noted 
above. Therefore, we treat star forma t ion in a way similai0 
to that introduced in INavakshin et all (|2007f ). Our approach 
is dif f erent from the industr y practice (e.g . . ISpringel et al.l 
l2005l ; iBooth fc Schavel 120091 ; iDubois et all |2012l) in moch 
elling star formation in cosmological simulation, where (i) 
a threshold density for star formation is introduced (rang- 
ing from 10 4 to 1 Hydrogen atom per cm 3 , typically); and 
(ii) once gas density exceeds the threshold stars are intro- 
duced stochasti cally following th e 3D form of the Schmidt- 
Kennicutt law (Kc nnicuttlll998r ). This method for treating 
star formation is reasonable on galaxy-wide scales, provided 
that the threshold density exceeds the mean gas density by 
a large factor. However, the method will fail in the inner 
parts of bulges where the tidal density, p t = o 2 /(2nGR 2 ) ~ 
W~ 22 R^p C g cm -3 ~ 100-Rj^ cm -3 is a strong function of 
radius R and may be higher than the threshold density. We 
also note that in the standard prescription, once the mean 
gas density exceeds the density threshold then stars are in- 
troduced everywhere in this prescription, even in gas regions 
that are not gravitationally self-bound, e.g., in the hot in- 
terclump medium. 

Our method avoids these potential pitfalls. Firstly, we 
require that the gas is self-bound before star formation is 
initiated. This is done by requiring that the Jeans mass, 
Mj ~ 2.6 x I0 e lf /2 n2 M Q , exceeds the minimum re- 
solvable mass of the simulation, i.e., that of 40 SPH par- 
ticles, m rcs = 2 x 10 5 M Q ; here T 4 = T/(10 4 K) and 
m = n/(10 4 cm -3 ). This criterion stops spurious star forma- 
tion in dense but hot gas (if such is formed in a strong shock, 
for example). Secondly, we require that the gas density 
should exceed the tidal density by a large factor, A = 200. 
The latter is a free parameter, but we tested that the re- 
sults change insignificantly as long as A is large, by running 
a set of identical tests in which A was varied from 2 to 200. 
Physically, this second condition prevents star formation in 
regions that appear to be self-gravitating but are not dense 
enough to withstand galactic or bulge tides. 

We are therefore confident that our methods correctly 
capture the onset of star formation, and do this better than 
the industry standard. However, due to our focus on the ef- 
fects of AGN feedback, we neglect star formation feedback 
within the gas clumps. The result of this approximation is 
that once star formation occurs within the clump it is likely 
to consume all the available gas, resulting in up to 100% ef- 
ficiency of star formation. This is most likely incorrect. Star 
formation feedback could disrupt the star forming clumps or 
heat them up, limiting star formation efficiency (although 
we note that the star forming phase of our simulations is 
quite short, from a few Myr to ~ 10 Myr, so only the most 
massive stars would have been able to produce supernovae 
quickly enough). Concluding this discussion, our methods 



1 Although due to a much poorer mass resolution here we do not 
attempt to model detail of sub-Sola r protostar e volution, such as 
potential mergers of "First Cores" llLarsorJll969l) or gas accretion 
onto the sink particles. Such microphysics of star formation is 
beyond the numerical realm of this paper. 
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Figure 1. The gas column density and velocity map for simula- 
tion M5, 11.8 Myr after the beginning of the simulation. Colour 
represents projected gas column density (scale at the bottom) 
and arrows indicate gas velocity. The arrows are logarithmically 
scaled to the longest one, so they only give qualitative information 
about gas motions. 



correctly predict when star formation occurs, but we prob- 
ably strongly over-estimate the amount of stars formed. We 
shall keep this conclusion in mind and come back to it in 
section 3. 

Here we present three simulations, Ml, M2 and M5, 
where the number indicates the mass of the SMBH in units 
of 10 s Mq. Therefore in the first simulation, we would ex- 
pect the shell to stall if it was driven predominantly by the 
momentum, rather than energy, input from the quasar. In 
the other two simulations, the shell should escape to infinity. 
We present the results below, drawing particular attention 
to the dynamics of the self-gravitating clumps formed in the 
three simulations. 



3 RESULTS 

Figure [T] shows the column density for the simulation M5 
and velocity vectors projected on to the plane of the figure 
for gas at 11.8 Myr after the quasar switches on. We see that 
while the bulk motion of the gas is radially outwards, there 
is a lot of structure forming as the shell fragments. This is es- 
sentially the same result as found by NZ12: the cooling time 
in the shocked ambient medium is t coo i ~4x 10 5 U8-Rk pc y r i 
where v$ = t> c /1000 km/s; this is a few times shorter than 
the flow timescale R v c ~ 10 6 yr, so the outer shell cools 
rapidly and is strongly compressed. In this simulation, the 
outflow does not stall, but rather keeps expanding and ac- 
celerating up to a radial velocity v r ~ 2000 km/s. This is 



because the SMBH mass is above the critical M a limit (eq. 
[1} and so even gas close to the SMBH is driven outside the 
cooling radius, where it can be accelerated further by the 
energy input. The network of filaments produces density en- 
hancments that become self-gravitating. 

In Figure [5] we show a narrow wedge projection (see 
NZ12 for more detail) of the simulation M2 at 7.09, 8.27 
and 9.45 Myr. White squares show star particles. Note that 
they form filaments pointing mainly in the radial direction, 
similar to the "fingers of creation" found in NZ12. Denser gas 
regions experience a lower outward push per unit mass, and 
therefore lag slighly behind lower density gas. These dense 
filaments eventually collapse into self-gravitating clumps, 
which are then observed to lag behind the shell. Clumps 
born earlier on are born with lower radial velocity since the 
gas shell is being accelerated outward due to its decreas- 
ing mean projected column density. Therefore dense gas re- 
gions collapsing over the course of the simulations result in 
"streamers" of star particles clearly visible in Figure 2. For 
the simulation M2 in particular, the first clumps move with 
mean radial velocities v» ~ 200 km/s, slightly lower than the 
bulge escape velocity v csc — 400 km/s (see also red points 
in Figure [3}- 

Figure [3] shows the radial velocity versus radius for star 
particles formed in all three simulations. The snapshots are 
taken at 16.5 Myr for Ml and at 9.45 Myr for the other two 
simulations. The dashed line in the simulation shows the 
escape velocity as a function of radius, assuming the effective 
radius for the model galaxy of 20 kpc. In simulation Ml, the 
shell fragments, the dense clumps stall and then fall back 
on to the SMBH. For this reason, there are both positive 
and negative radial velocities for stars in this simulation. 
Almost all of the stars formed in this simulation (blue dots) 
are bound to the host galaxy. 

For the more massive SMBHs, the shell is continuously 
driven out. As the shell is of a finite mass, its mean pro- 
jected column density decreases as it expands, and so quasar 
outflow can accelerate it to velocities larger than would 
be c alculated analyti cally for a shell of infinite radius (see, 
e.g.. lKing et al.ll201ll . for derivation). As already explained 
above, the collapsing star-forming clumps do not experience 
the same outward acceleration from the quasar outflow and 
decouple from the shell. Star particles born earlier are born 
with a lower outward velocity, therefore they lag behind 
those that are born at later time. As a result the dispersion 
in clump velocities increases with time, which is already ev- 
ident from Figure [2] Figure [3] shows that some of the stars 
formed in simulation M2 are bound to the galaxy (those be- 
low the dashed curve); they move on very elongated, nearly 
radial orbits. Stars formed later have large enough radial ve- 
locities to escape the potential completely. In simulation M5 
the whole of the shell is accelerated to velocities above the 
escape velocity before the stars are born; all of the newly 
made stars are on escaping trajectories. 

A large fraction, 30 — 70%, of the initial mass of the 
shell is converted into star particles in each of the three 
simulations. As discussed in section such a large star for- 
mation efficiency is clearly an overestimate. Star formation 
feedback, neglected here, is likely to reduce the total mass 
of the newly made stars. Nevertheless, even if star forma- 
tion efficiency drops by an order of magnitude, e.g., to 5%, 
the total mass of stars formed on the radial outward tra- 
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Figure 2. A wedge-slice through simulation M2 at times 7.09, 8.27 and 9.45 Myr (left to right panels, respectively). Colours show 
projected gas surface density, white squares represent star particles. Stars form from overdensities in the fragmenting shell, but their 
clusters are stretched out because the stars decelerate while the shell moves outwards with approximately constant velocity. Arrows are 
gas velocity vectors, giving qualitative information about the turbulence in gas motion induced by quasar feedback. 
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Figure 3. Radial velocity of stars formed in simulations Ml, 
M2 and M5 versus their instanteneous radial coordinate, R. The 
dashed line shows the escape velocity as a function of R. The star 
forming shell created by the lightest black hole stalls at a radius 
of R ~ 0.3 kpc, so all of the stars formed in this simulation are 
bound to the host. 



4 DISCUSSION 

The simulations presented here show that, in the setup ex- 
plored, quasar feedback can expel gas from the bulge of its 
host galaxy with high velocity and simultaneously compress 
it to high densities. These two processes result in forma- 
tion of hyper- velocity stars (HVS). The total mass of these 
stars is unclear at the moment, since our simulations neglect 
star formation feedback (see [}2]and[3]). In addition, here we 
explored the simplest possible initial conditions - a spher- 
ically symmetric gas distribution, with the density profile 
p oc 1 / R 2 . Simulations of AGN feedback with more realistic 
non-spherical initial conditions show that the quasar outflow 
driven bubble may burst in directions of the lowest column 
depth of the ambie nt gas, e.g., perpendicular to the planet 
of the host galaxy l|Navakshin fc Z ubovas 2012). Feedback 
in such a geometry may be ( Zubovas fc Nayakshin! I2012T ) 
the reason behind formatio n of the "Fermi Bubbles" in the 
Milky Way (|Su et al.ll2oiol ). In such geometries the fraction 
of ambient gas put on hyper-velocity orbits is smaller, and 
therefore we expect that the total mass of the HVS would 
be also smaller. 

Therefore, while we are confident that AGN feedback 
can launch HVS, we are currently unable to quantify how 
often this occurs in more realistic situations, and how large 
the fraction of gas turned into stars is. With this caviat in 
mind, we now discuss possible observational implications of 
HVS made by AGN feedback. 



jectories would be between a few times 10 M© to 10 Mq . 
We also note that star formation efficiency is expe cted to 
be higher in compressed systems dKeto et al.l 120051 ). Thus, 
provisionally, we may expect at least ~ 10 8 M© of stars to 
be formed. 



4.1 Radial mixing of stars and clouds 

Stars formed in the shock driven outward by the quasar with 
velocities below the escape velocity, such as those shown 
with the red symbols that are below the dashed curve in 
Figure [3] are destined to fall back to the bulge. These stars 
are on very elongated nearly radial trajectories, and spend 
most of the time at large radii, i.e., R ^ a few to a few 
tens of kpc. These stars are thus unusual in that they were 
born from the gas initially filling the inner region of the 
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galaxy, the bulge, but due to their large radial velocity kicks 
now belong to the galactic halo. Such halo stars could thus 
stand out due to their highly radial orbits and (likely) higher 
metallicity compared to a typical low metallicity halo star. 
The most extreme of these stars, in the outskirts of the dark 
matter halos of their host galaxies, may provide a significant 
contribution to the intergalactic infrared background. It was 
recently shown that a 0.1% fraction of total galaxy light de- 
posited into the outskir ts of halos could exp lain most of the 
background radiation (jCoorav et al.l l2012h : while we can- 
not quantify the total luminosity of stars expelled to such 
distances in our model, our results are consistent with this 
requirement. 

The gas clumps expelled from the bulge do not nec- 
essarily all form stars. Some of these clumps may remain 
marginally bound, and collapse on a timescale longer than 
the cloud ejection timescale. In that case, the gas clumps can 
mix with the halo gas, joining the high velocity cloud (HVC) 
population. These clouds are observed all around the Milky 
Way, have high peculiar veloc ities and metallicities similar to 
or slightly low er than Solar l|Wakker fc van Woerden|[l997l ; 
IWakkeij[200ll ). Their evolutionary timescales are estimated 
at a few tens to several hundred Myr, suggesting continu- 
ous replenishment, su ch as condensation fro m tidal streams 
JPutman et al.ll2p03bl ) or a galactic fountain l|Bregmanlll980l ; 
iKwak et al.ll2009r ). Our simulations suggest that quasar out- 
flow, long extinguished by now, could also contribute to pro- 
duction of HVCs. 



4.2 Bulge disruption 

It is usually assumed that gas is a small fraction of the total 
mass of the host, with the rest consisting of dark matter 
and stars. This is clearly correct for present day galaxies, and 
must be globally true for galaxies at birth since the universal 
baryon fraction is only / g ~ 0.16. However, gas can collapse 
to the centre of the potential if its angular momentum is 
sufficiently low and so dominate the potential in the central 
region of the host. In this case, if most of the gas is driven 
rapidly out of the host, the remaining mass may be insuffient 
to bind the bulge gravitationally. It then has to expand and 
perhaps blend with the rest of the galaxy, adding stars to 
the halo, with faster stars leaving the host galaxy altogether. 

In the extreme cases this should result in an SMBH 
without a clear parent bulge or an SMBH that appears to be 
far more massive than it should be based on the A/bh — -A/b u i go 
relation. This cannot be a common outcome: based on Fig- 
ure 3, we see that the SMBH mass needs to be somewhat 
above the M a value to drive the cold shell out before the lat- 
ter collapses into stars. This happens even in a spherically- 
symmetric geometry, which is more conducive to large out- 
flows; NZ12 found that for more realistic non-spherical ge- 
ometries, AGN feedback in the gas rich epoch tends to com- 
press and trigger star bursts in cold gas without driving the 
shell away. 

Nevertheless, bulge destruction or erosion (partial de- 
struction) by the SMBH feedback expelling a large fraction 
of the bulge mass may occasionally happen, especially for 
active galaxies in the early Universe and/or violently merg- 
ing systems. This may be relevant for the re cently observed 
popu l ation of SMBHs in bu l geless galaxies (fSimmons et al.l 
120121 ; iMcAlpine et all l201ll ; lArava Salvo et all l2012h and 



well as several SMBHs with masses significantly above 
those predicted by t he black hole - bulge mass correlation 
i|Bogdan et al.ll2012l ). 



4.3 Type II supernovae on galactic outskirts 

The distance travelled by a star moving with radial velocity 
v — 1000i>8 km s _1 for time t — 10 7 t7 yrs is 



AR = vt « 10 kpc v 8 t 7 



(2) 



Stars with Solar or super-Solar metallicity that are more 
massive than M ~ 8Mq end their main seque nce life in a 
Type II supernova explosion (He ger et al.l 120031 ). The main 
sequence lifetime of stars of mass M — 8 Mq is approxi- 
mately 100 Myr, while the most massive stars, M ^ 50 Mq 
live only ~ 4 Myr. Therefore we can expect that mas- 
sive stars launched by the quasar outflow produce Type 
II supernovae at a distance from a few to few tens of kpc 
away from the galaxy's center. Additionally, long duration 
GRBs are expected from this flung-out population of stars 
l|Wooslev fc BloomlbOOq ). potentially explaining the GRBs 
that are known to be offset, some by a s much as > 10 kp c, 
from the centres of their host galaxies (|Bloom et al.| [2002) . 



4.4 Reionisation and metal enrichment of the 
Universe 



iBouwens et~aH il201ll ) show that stars formed in galaxies at 
around the reionization epoch can provide sufficient amount 
of ionizing photons provided that they escape into inter- 
galactic space efficiently enough. In particular, these authors 
estimate that an escape fraction of 20% is required for con- 
sistency with the WMAP results at the 2a level. This is 
much higher than the esca pe fraction of only ~ 1 — 2% as 
found for the Milky Way l|Putman et all l2003aT ). The dis- 
agreement is all the more striking because the present day 
Milky Way contains very little gas per unit total mass of 
the Galaxy compared with the early star forming galaxies, 
which were also probably much more compact given that ob- 
servations show that gala xies grew inside out at le ast from 
z ~ 2 to the present day l|van Dokkum et al.ll2010T l. 

It is therefore interesting to note that the massive stars 
launched at high velocities, v ^ v asc , find themselves out- 
side the cold high density gas environment of the inner 
galaxy. Ionizing photons emitted by these stars have orders 
of magnitude smaller average gas column depth to navigate 
through to escape from the galaxy's halo. Note that while 
massive stars are trailing the cold shell that created them 
earlier on, the shell is a web of thin dense filaments with most 
of the sky covered by a much lower density material (see Fig. 
[T]and also Fig. 2 in NZ12). Therefore, the escape fraction 
of ionising radiation for these stars may be much closer to 
unity than for the stars residing in the inner galaxy. The 
Hyper Velocity Stars (HVS) formed during quasar outflow 
ejection may thus be far more efficient in terms of reioniza- 
tion than the "normal" galactic stars, although it remains 
to be seen if a sufficient mass of HVS can be created to 
contribute significantly to the reionization of the Universe. 

Lower mass stars ejected of their host galaxies at high 
velocities may be contributing to metal enrichment of the 
gas outside of galaxies. Their longer main sequence lifetime, 
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t ~ 10 9 yrs, implies that they may travel as far as ~ 1 Mpc 
(see eq. [2]) from their hosts before ejecting their envelopes. 
Recent observations have discovered low-mass stars at dis- 
tances of up to 2 Mp c from the centre of the Milky Way 
|Palladino et al.ll2012T l. These stars might have been ejected 
by three-body interactions between a binary star and th e 
central SMBH of the Milky Way fe.g jYu fe Tremainei20o3 ). 
but without knowledge of their orbits, it is impossible to con- 
firm or disprove this scenario. The existence of these stars is 
also consistent with our model; they might have been ejected 
by a quasar outflow from the Milky Way several Gyr ago. 
Importantly, this "stellar outflow" may actually reach fur- 
ther into the intergalactic space than quasar gas outflows. 
Here we found that stars are lagging behind gas outflows 
while quasar outflow provides an outward acceleration for 
the gas shell; but eventually the quasar must turn off, and 
the gas outflow must decelerate faster than a HVS because 
the outflow interacts with intergalactic medium by a shock 
whereas the stars continue to travel ballistically. 



5 SUMMARY AND CONCLUSION 

Using numerical simulations of SMBH feedback upon gas 
in its host galaxy, we have shown that quasars may not 
on ly trigger star formation in their host bulges (as shown 
bv lNavakshin fe Zubovas]|2012h . but also expel a fraction of 
these newly-formed stars from the host galaxy. Expulsion 
requires a large SMBH mass (M > M a ), so may not be a 
common phenomenon. Further modelling is needed to de- 
termine the mass fraction produced in this HVS population, 
since its potential implications are wide ranging. 

Stars originating in the bulge launched at high veloci- 
ties on nearly radial trajectories may completely escape the 
galaxy; slower stars add a metal-rich subpopulation in the 
region and may contribute to the intergalactic light. High 
metallicity HVS produce ionising radiation and so may con- 
tribute to the reionisation of the Universe; they also explode 
as supernovae, potentially contributing to Type II super- 
novae on galactic outskirts. Lower mass stars travel further 
out, enriching the extragalactic medium with metals. Fi- 
nally, a significant fraction of the galaxy bulge mass may 
become unbound due to such "stellar outflows", leading to 
SMBHs in bulgeless galaxies or SMBHs with masses much 
higher than expected from the black hole - bulge mass rela- 
tion. We conclude that quasar outflows shape the properties 
of their galaxies and surrounding medium in many ways, 
not just by removing gas and quenching star formation, as 
usually acknowledged. 
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